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ABSTRACT 

We observed the unidentified TeV gamma-ray source HESS J1804-216 with the CANGAROO-III 
atmospheric Cerenkov telescopes from May to July in 2006. We detected very high energy gamma 
rays above 600 GeV at the IOct level in an effective exposure of 76 hr. We obtained a differential 
flux of (5.0 ± 1.5stat ± l-6sys) X 10-i2(E/l TeV)-" cm-^s-^TeV"! with a photon index a of 2.69 ± 
0.30stat±0.34sys, which is consistent with that of the H.E.S.S. observation in 2004. We also confirm the 
extended morphology of the source. By combining our result with multi-wavelength observations, we 
discuss the possible counterparts of HESS J1804— 216 and the radiation mechanism based on leptonic 
and hadronic processes for a supernova remnant and a pulsar wind nebula. 

Subject headings: gamma rays: observations — ISM: individual(HESS J1804— 216, G8.7— 0.1) — 
pulsars: individual(PSR B1800-21) — X-rays: individual(Suzaku J1804-2142, 
Suzaku J1804-2140) 



1. INTRODUCTION 

A Galactic plane survey was performed in 2004 by 
the H.E.S.S. imaging atmospheric Cerenkov telescope 
(lACT) with a fiux se nsitivity of 0.02 c rab fo r gamm a 
rays above 200 GeV (jAharonian et al.l l2005aL l2006ar ). 
Fourteen new gamma-ray sources were detected at signif- 
icance levels above 4(T, and 11 of the sources either have 
no counterpart or possible counterparts with significant 
positional offsets. HESS J1804— 216 is one of the bright- 
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est, and its spectrum is softest in this survey; the fiux 
is about 0.25 crab above 200 GeV with a photon index 
of 2.72 ± 0.06. In addition, with a size of ~ 22 arcmin, 
it is one of the most extended TeV gamma-ray sources. 
The H.E.S.S. collaboration proposed two possible coun- 
terparts: the supernova remnant (SNR) G8.7— 0.1 and 
the young Vela-like pulsar B1800— 21. However, the TeV 
gamma-ray source does not coincide exactly with either 
of these. G8.7— 0.1 appears as a larger circular region 
with a diameter of 50 arcmin, and the geometric center 
(|Kassim fc Weilei1[r9 90bl has a large offset, ~ 11 arcmin, 
from the centroid of HESS J1804-216 to the northeast, 
while PSR B1800-21 has an - 11 arcmin offset to the 
west. 

After the detection of HESS J1804-216, an SNR, 
G8. 31— 0.09, was discovered at radio wavelengths 
(jBrogan et al.l 12006) to be located within the error cir- 
cle of HESS J1804-216, but with a smaller size, 5' x 
4'. Suzaku deep observations discovered two new X- 
ray sources, Suzaku J1804— 2142 (hereafter Srcl) and 
Suzaku J1804— 2140 (hereafter Src2), that are near the 
center of HESS J1804-216 ( Bamba et al.ll2007l) Chan- 
dra also detected these sources (|Kargaltsev et al.ll2007b[ ) . 
SWIFT found thre e faint X-ray sources in the region of 
HESS J1804-216 (iLandi et all 12001) . One of them is 
positionally coincident with a bright star, and another 
could also be associated with a star close to the bound- 
ary of the XRT error circle. The other positionally co- 
incides with Suzaku Src2. So, there remain 5 possi- 
ble counterparts: SNR G8. 7-0.1, PSR B1800-21, SNR 
G8.31-0.09, Suzaku Srcl, and Suzaku Src2. We briefiy 
describe these sources in the following paragraphs. 

SNR G8. 7-0.1: G8.7-0.1 is associated with the W30 
complex, which comprises extended radio emission with 
a number of superposed smaller discrete emission regions 
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(lAltenhoff et al.l 119781: iReich et all 11984 iHanda et all 
Il987( ). Radio recombination-line observations have been 
used to identify discrete sources as HII regions, and 
CO observations also show molecu lar gas to be asso- 
ciated with W30 (Blit z et al.lll982l ). lOjeda-Mav et all 
(j2002i) reported that massive star form ation may be oc - 
curr mg m molecular c louds i n W30. lOdegardI (|1986f ) 
and iKassim fc Weileil (|1990b( ) clearly established that 
G8.7— 0.1 was an SNR by the detection of a nonthermal 
extended radio emission. A i? 05'^ T" observation revealed 
diffuse X-ray emission only from the northern half of 
the r emnant in the 0. 1-2.4 keV band (|FinleY fc Qgelmanl 
[1991 . The distance to G8. 7-0.1 was estimated us- 
ing several methods. Based on kinematical distances 
to the HII regions associated with the SN R, the dis- 
tance was estimated to be 6 ± 1 kpc (Kassim fc Weileil 
ll990bD . iFinlev fc Qgelmanl (|1994D pointed out that more 
recent galactic rotation models applied to the H II regions 
suggest a near kinematical distance of about 4.8 kpc. 
They a lso estimated the distance based o n a Sedov so- 
lution (|Sedovlll959l: iHamilton et~al]|1983D from the ob- 
served X-ray temperature and the angular radius to de- 
rive 3.2 kpc < d < 4.3 kpc for an assumed initial energy 
of 10^^ ergs. In this paper, we adopt d = 4.8 kpc. They 
also estimated the age of the SNR, based on a Sedov 
solution from the X-ray observation, to be 1.5-2.8 xlC* 
years under the assumption of an initial energy of 10^^ 
ergs, and 2.7-3.9x10^ years under the assumption of 
a distance of 6 kpc. Similarly, we estimated an age of 
2.2-3.1 xlO^ year s unde r the assumption of a distance of 
4.8 kpc. lOdegardI (|1986l ) represented an age of 1.5 x 10^ 
years from the relation between the age and the surface 
brightness in the radio band. 

PSRB1800-21: The young Vel a-like pulsar B1800-2 1 
was found in a radio observation (jClifton fc Lvnel[l986f ). 
The offset from the centroid of HESS J1804-216 is 
~11 arcmin. The spin period, the spin period deriva- 
tive, and the characteristic age are P = 133.6 ms, 
P = 1.34 X IQ-i^s/s, and r^. = P/2P = 15.8kyears, 
respectively (|Brisken et al.l [20061 ). The resulting spin- 
down luminosity is E = in'^IP/P^ = 2.2 x 10^^{I/W^^) 
ergss"^, where / is the moment of inertia in units 
of gcm^. The distance was estimated to be 3.9 kpc 
from the pulsar's dispe rsion measure of 233.99 pccm"'^ 
(jTavler fc Corded I1993D. The n ewer Cordes & Lazio 
NE2001 model (^Cord es fc Laziol |2002^ gives a similar 
distance of 3.84i[]j^ kpc. We adopt d = 3.84 kpc 
throughout this paper. The association between PSR 
B1800— 21 a nd G8.7— 0.1 has been d i scussed in several 
papers (e.g. iFin ev fc Qgelmanl [l99l iFrail et al.l Il994l : 
IKassim fc Weileil Il990air~ ^However, a re cent proper- 
motion measurement (jBris kcn et al.] |2006| ) has shown 
that the pulsar was born outside the currently observed 
SNR, and that is moving more nearly toward the cen- 
ter of G8.7— 0.1, rather than away from it, which makes 
their association very unlikely. Based on a 10 ks ob- 
servation with the ROSAT PSPC, Finley fc Ogelman 
(1994) reported a faint X-ray source near the radio 
pulsar position, and at tributed this emi ss ion to PSR 
B1800-21. Rec ently iKargaltsev et all (|2007al ) and 
ICui fc Konopelkol (|2006l ) reported that an X-ray nebula 
arou nd the pulsar was detected with Chandra. Addition- 
ally, IKargaltsev et al.l ((20073) reported that the X-ray 



nebula has two structures: a brighter compact 7" x4") 
component (the inner pulsar wind nebula (PWN)) and an 
extended (~ 12") fainter emission component (the outer 
PWN). These are asymmetric to the pulsar position and 
extended toward HESS J1804-216. 

SNR G8.31-0.09: SNR G8.31-0.09 was found in a 
90 cm multi-configuration Very Large Array survey of the 
Galactic plane. The size is 5' x 4', and the offset from 
the centroid of HESS J1804-216 is 7 arcmin. The mor- 
phology is shell-like and the spectral index is a^. — —0.6 
for cx I/"- (Brogan et al. 2006). 

Suzaku Srcl: The offs et from the c e ntroid of HESS 
J1804-216 is 3 arcmin. iBamba et aH ()2007l ) reported 
that Srcl is point-like or compact compared to the spatial 
resolution of Suzaku with a half-power diameter of about 
2 arcmin. However, Kargaltsev et al. (2007b) reported 
that Srcl is extended or multiple (1.5' — 2') with a Chan- 
dra observation. The Suzaku spectrum was fitted with 
an absorbed power-law model ()Bamba et all [2007). The 
best-fit absorbing column is consistent with the Galactic 
hydrogen column in that di rection. Since the ph oton in- 
dex of —0.3 ±0.5 is very flat. lBamba et all (|2007| ) suggest 
that this source is likely to be a high-r nass X-ray binary 
(HMXB). In the Chandra observation ([Kargaltsev et al] 
l2007br ). no spectral fitting was able to be performed be- 
cause of the low signal-to-noise ratio. The unabsorbed 
flux, which was estimated from the Chand ra observation 
using the best-fit parameters reported by iBamba et al.l 
(|2007h . is a factor of « 1.7 smaller than that reported 
with the Suzaku observation ofl Bamba et all ([20071 ) . The 
difference could be due to unaccounted systematic er- 
rors, or the variability of the source, which sup ports the 
HMXB interpretation ([Kargaltsev et alJl2007H ). 

Suzaku Src2: The offset from t he centroid of HESS 
J1804-216 is 1.8 arcmin. Though [Bamba et all ([2007f ) 
reported that Src2 is extended or mult i ple, ba sed on the 
Suzaku observation, IKargaltsev et al.l ( 2007bf l reported 
that Src2 is point-like in a Chandra observation. This 
could mean that the more sensitive Suzaku observations 
have detected a fainter extended PWN component. The 
Suzaku spe ctrum was fitted w ith an absorbed power- 
law model ([Bamba et al.|[2007[ ). and the best-fit absorb- 
ing column is about an order-of-magnitude higher than 
the expected Galactic column. This implies that Src2 
is embedded in dense gas. The Chandra spectrum was 
also well-fitted with an absorbed power-law model, and 
the obtained absorbing column densi ty is a factor of 
2-3 la rg er than the Gala c tic colu mn ([Kargaltsev et al] 
I2007bl) . IKargaltsev et all ([2007bO stated that the large 
absorption suggests that Src2 is located within (or even 
beyond) the Galactic Bulge, or it shows an intrinsic 
absorption, which is often se e n in th e X-ray spectra 
of HMXBs. IKargaltsev et all ([2007bO also reported a 
marginal pulsation of 106 s in Src2, whic h supports an 
HMXB interpretation. On the other hand, IBamba et alJ 
([2007) suggested that Src2 is a PWN or a shell-like 
SNR because of the extended morphology observed with 
Suzaku and the best-fit photon index of 1.7 (0.7-3.1). 

None of the above five sources morphologically match 
HESS J1804-216. The counterpart is therefore still un- 
known. In this paper, we present TeV gamma-ray ob- 
servations of HESS J1804-216 with the CANGAROO- 
III telescopes and discuss the radiation mechanism and 
the counterpart by combining our result with multi- 
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wavelength observations. 

2. CANGAROO-III OBSERVATIONS 

CANGAROO-III is an array of four lACTs, located 
near Woomera, South Austraha (136°47'E, 31°06'S, 
160 m a.s.l.). Each telescope has a 10 m diameter re- 
flector made up of 114 segmente d FRP spherical mir - 
rors mounted on a parabolic frame (|Kawachi et al.|[200lD . 
The telescopes are situ ated at the corners of a diamond 
with sides of ~ 100 m ()Enomoto et al.ll2002D . The old- 
est telescope, Tl, which was the CANGAROO-II tele- 
scope, was not used due to its smaller FOV and higher 
energy threshold. The imaging camera systems on the 
other three telescopes (T2, T3 an d T4) are identica l, 
with 427 PMTs and a FOV of 4.0° (|Kabuki et al.l[200l ). 
The PMT signa ls were recorded by charge ADCs and 
muhi-hit TDCs (jKubo et al.ll2001[ ). The observations of 
HESS J1804-216 were made from May to July, 2006, us- 
ing the 'wobble' mode in which the pointing position of 
each telescope was shifted in declination by ±0.5° from 
the centroid of HESS J1804-216. The mean zenith an- 
gle of the observation was 21°, and the total observation 
time was 86.8 hr. We used the 3-fold coincidence data 
taken at zenith angles of less than 40°. To trigger data 
recording, an individual telescope was required to have 
more than four pixels receiving over 7.6 photoelectrons 
within 100 ns (local trigger), with a global trigger system 
to deter mine the coincidence o f any two of the three tele- 
scopes (jNishiiima et al.l[2005H . We rejected data taken 
under bad weather conditions in which the shower event 
rate was less than 6 Hz. Taking into account the DAQ 
dead-time, the effective live time was calculated to be 76 
hr. 

3. DATA REDUCTION AND ANALYSIS 

The basic analysis pro cedures are described in detail 
in lEnomoto et all (|2006a[ ) and lKabuki et all (|2007D . Us- 
ing calibration data taken daily with LEDs, the recorded 
charges of each pixel in the camera were converted to the 
number of photoelectrons. At this step we found 7 bad 
pixels out of 427 pixels for T2, 5 for T3, and 1 for T4, 
due to their higher or lower ADC conversion factors in 
these observations. These bad pixels were removed from 
this analysis, which was also reflected in the Monte Carlo 
simulations. After that, every shower image was cleaned 
through the following CANGAROO-III standard crite- 
ria. Only pixels that received >5.0 photoelectrons were 
used as "hit pixels". Then, five or more adjacent hit pix- 
els, with arrival times of within 30 ns from the average 
hit time of all pixels, were recognized as a shower cluster. 

Before c alculating im age moments — the "Hillas pa- 
rameters" dffillas' IQSSIL — we app lied the "edge cut" 
described in^nomoto e t al.l (|2006bl ). We rejected events 
with any hits in the outer-most layer of the camera. The 
orientation angles were determined by minimizing the 
sum of the squared widths with a constraint given by 
the distance predicted by Monte Carlo simulations. 

We then a p plied the Fisher Disc riminant method 
(iFished Il936t lEnomoto et~al] l2006aD with a multi- 
parameter set of P = (W2, W3, W4, L2, ^3, ^4), where 
W and L are the energy corrected width and length, and 
the suffixes represent the telescope IDs. The Fisher Dis- 
criminant (FD) is defined as FD = a ■ P, where a is 
a set of coefficients mathematically determined in order 



to maximize the separation between two FDs for gamma 
rays and hadrons. 

For a background study we selected a ring region 
around the target, 0.3 < 9^ < 0.5 deg^, where 9 is the 
angular distance to the centroid of HESS J1804— 216 re - 
ported by the H.E.S.S. group tAharoni an et al.ll2006af ). 
and obtained the FD distributions for the background, 
Fb, and Monte Carlo gamma rays, Fg. Finally, we could 
fit the FD distributions of the events from the target 
with a liner combination of these two components. The 
observed FD distributions, F, should be represented as 
F — aFg + (1 — a)Fi,, where a is the ratio of gamma-ray 
events to the total number of events. Here, only a was 
optimized. This analysis method was verified by an anal- 
ysis of the Crab nebula data taken in December, 2005. 

The reflectivities of each telescope, which were used 
in the Monte Carlo simulations, were monitored every 
month by a muon ring analysis of a calibration run taken 
individually by each telescope. We obtained relative 
light-collecting efficiencies with respect to the original 
mirror production times of 0.60, 0.60 and 0.65 for T2, T3 
and T4, respectively. Throughout this analysis, we used 
the Monte Carlo simulations for gamma rays assuming a 
point-source. 

4. RESULTS 

The obtained 9^ plot is shown in Fig. [1] with the point 
spread function (PSF) of our telescopes, 0.23° (68% con- 
tainment radius). The numbers of excess events that we 
detected above 600 GeV were 512 ± 61 within 9'^ < 0.06 
deg^, based on the assumption that it was a point source, 
977 ± 94 within 9"^ < 0.17 deg^, which corresponds to 
that used in the spectral analysis by H.E.S.S. (and tak- 
ing into account the difference between our PSF and 
that of H.E.S.S.), and 1389 ± 126 within 9"^ < 0.3 deg^. 
The TeV gamma-ray emission is extended, and the mor- 
phology of gamma-ray-like events, derived with boxcar 
smoothing with each pixel replaced by the average of its 
square neighborhood, is shown in Fig.[2l The number 
of excess events was individually estimated by the FD- 
fitting method in each 0.2° x 0.2° sky bin. When we 
evaluated the outer regions {9^ > 0.6 deg^), we had to 
consider gradual deformations of the FD distributions at 
larger angular distances from the target. Therefore, we 
selected an annulus with radii 0.2° < r < 0.4° centered 
on the evaluated region as the background. For the inner 
regions, 9^ < 0.6 deg^, the events in 0.3 < 9"^ < 0.5 deg^, 
excluding the evaluated region, were adopted as a back- 
ground. The intrinsic extent of the TeV gamma-ray 
emission was estimated by a 2D Gaussian fit on our un- 
smoothed excess map. The intrinsic deviations along 
the Right Ascension and Declination axes were calcu- 
lated to be 0.160° ± 0.005° and 0.274° ± 0.011°, respec- 
tively. The best-fit centroid position was obtained (R.A, 
dec [J2000])=(271.079°, -21.727°). Th e off'set from the 
best-fi t position reported by H.E.S.S. (jAharonian et all 
[2006al) is (AR.A, Adec)=(-0.053° ± 0.007°, -0.026° ± 
0.013°). The offset is not significant given our PSF. 

Figure [3] represents a reconstructed VHE gamma- 
ray differential spectrum compatible with a single 
power-law: (5.0 ± 1.5stat ± l-Gsys) x 10-i2(e/1 TeV)"" 
cm-^s-^TeV-i with a photon index a of 2.69±0.30stat± 
0.34sys. To obtain the spectrum, we used a cut of 
9^ < 0.17 deg'^. The relevant systematic errors are due to 
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Fig. 1. — Number of excess events as a function of the 
squared angular distance. Here, 0° corresponds t o the centroid of 
HESS J1804-216 reported by the H.E.S.S. group l lAharonian et al.l 
[2006a). The squares show the CANGAROO-HI data points. The 
circles show the normalized H.E.S.S. data points. The hatched 
histogram represents our PSF for a comparison. 

the atmospheric transparency, night sky background fluc- 
tuations, uniformity of camera pixels, and light-coUecting 
efficiencies. In addition, the signal integrating region 
was changed from 9"^ < 0.17 deg^ to 0.3 deg^, and the 
difference in fluxes was incorporated in the systematic 
errors. The TeV gamma-ray extension and the flux ob- 
tained by CANGAROO-III were consistent with those by 
H.E.S.S. Our result indicates that the TeV gamma-ray 
emission was unchanged between the H.E.S.S. observa- 
tions in 2004 and ours in 2006. 

5. DISCUSSION 

As described in § 1, there are 5 possible counterparts: 
SNR G8.7-0.1, PSR B1800-21, SNR G8.31-0.09, 
Suzaku Srcl, and Suzaku Src2. Figures 2] and [H] show 
the spectral energy distribution (SED) of these counter- 
parts and HESS J1804— 216 including the results of this 
work. We extend the introductions to these sources given 
in § 1 with a more detailed summary of the sources and 
t heir characteri s tics in the following paragraphs. 

iBamba et al.l ()2007[ ) state that Suzaku Srcl and Src2 
are physically associated with HESS J1804— 216. They 
suggested that X-rays and TeV gamma-ray emission 
could co me from an SNR shock, based on a model pro- 
posed bv lYamazaki et ahl ()2006D : in an old SNR with an 
age of ~ 10^ years, primary electrons have already lost 
most of their energy, and only nucleonic cosmic rays re- 
main. Additionally, the old SNR shock colliding with a 
giant molecular cloud (CMC) can emit hard nonthermal 
X-rays from secondary electrons and strong TeV gamma 
rays from shock accelerated protons through 7r° decay. 
This model can explain the observed large ratio of the 
TeV gamma-ray to X-ray flux (factors of ~100). The 
large X-ray absorption of Src2 also supports this sce- 
n ario. 

iKargaltsev et al.l (|2007bD suggested that a PWN is 
the source of HESS J1804-216, like HESS J1825-137 
or Vela X. HESS J1825-137 is hkely to be associated 
with the PWN G18.0-0.7 around the Vela-like pulsar 
B1823— 13. The TeV gamma-ray emission detected with 
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Fig. 2. — Smoothed morphology of gamma-ray-like events with 
our PSF of 0.23° radius. Dashed contours show the VHE gamma- 
ray emission seen by H.E.S.S. ( Aharonian et al. 2006a). The thick 
solid contours (green) show the 20 cm radio emission from G8.7— 0.1 
recorded by the VLA (White et al. 20051). The thin solid contours 
(white) show the X-ray emission detected by the ROSAT satel- 
lite l IFinlev fc Ogelm an||199 4|). The solid circle indicates the posi- 
tion orGSTsT^aOgTBroganeFn^ '2^^ The cross indicates the 
PSR B1800-21 position (Brisken et al. 2006). The triangle and 
the square indicate the position of Suzaku Srcl and Suzaku Src2, 
respectively. 
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Fig. 3. — Differential flux. The squares and circles show the 
CANGAROO-III and the H.E.S.S. data points, respectively. The 
best-fit power-law is also shown by the solid and dashed line from 
this work and from H.E.S.S., respectively. 

H.E.S.S. covers a much larger area than the X-ray emis- 
sion from G18. 0—0.7, extending up to 1° southward from 
the pulsar ijAharonian et al.]|2005bD . However, both the 
TeV gamma-ray and the low surface-brightness X-ray 
emission have similarly asymmetric shapes, and they are 
offset in the same direction with respect to the pul- 
sar position. A similar picture is ob served around the 
Vela pulsar (| Aharonian et al.) I2006bf ). These phenom- 
ena can be exp lained by the "crushed PWN" hypothesis 
(|Blondin et al.ir2001 1: on a time scale of ~ 10"' years, the 
reverse SNR shock front propagates toward the center of 
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the remnant, where it crushes the PWN, and asymme- 
tries in the surrounding interstellar medium give rise to 
an asymmetric shape and offset of the PWN relative to 
t he pulsar and exp lo sion sit e. 

Kargalt sev et al.l ()2007bl ) considered the possibility 
that Suzaku Srcl or Src2 are a PWN powering HESS 
J1804— 216. However, the 3.24 s time resolution of the 
Chandra ACIS observation precludes a search for the 
subsecond pulse periods expected for a young pulsar. 
Therefore, there is no strong evidence to support it at 
this point. They also suggested that PSR B1800-21 is 
associated with HESS J1804-216. Its asymmetric PWN 
component e xtending toward HESS J 1804— 216, detected 
by Chandra ([Kargaltsev et al]|2007af ). shows a hint of the 
association, but the sensitivity of the Chandra observa- 
tion was possibly insufficient to detect the PWN beyond 
15"-20" from the pulsar. Additionally, they pointed out 
that the extended morphology of HESS J1804-216 ar- 
gues against the HMXB interpretation because of the 
weak observational evidence for HMXBs producing ex- 
t ended TeV gamma-ra y emission. 

iFatuzzo et all (|2006D concluded that PSR B1800-21 
cannot account for the spectrum of HESS J1804— 216, 
and G8.7— 0.1 is probably the source of the TeV gamma 
rays. However, they considered only a pion-decay model 
without any consideration of the inverse Compton pro- 
cess, and for PSR B1800— 21 they only considered the 
acceleration of charged particles across voltage drops in 
the relativistic winds near the light cylinder. Addition- 
ally, they did not take into account the possibility of 
other sources besides PSR B1800-21 and G8. 7-0.1. 

Based on the above discussions and the SED, we now 
discuss the radiation mechanism of HES S J1804— 216 
and i ts counterpart. For Suzaku Srcl, iBamba et al.l 
(|2007f) reported a very hard photon index of —0.3 ± 0.5, 
which corresponds to the power-law index of electrons 
of —1.6 ± 1.0. This unusual value indicates that the X- 
ray emission is not synchrotron rad iation. Theref o re, we 
do not discuss the SED of Srcl. iBamba et all (|2007D 
suggested that Srcl might be an HMXB. If it were, the 
association between Srcl and HESS J1804— 216 wo uld be 
unlikely as discussed in Kargaltsev et all (|2007bD . The 
results of independent observations of the H.E.S.S. and 
CANGAROO-III telescopes show that HESS J1804-216 
is quite extended (~ 0.4°). This precludes the possibil- 
ity that HESS J1804-216 is an Active Galactic Nucleus. 
The plausible candidates seem to be an SNR or a PWN. 

5.1. SNR origin 

The SNR candidates associated with HESS J1804-216 
are SNR G8. 7-0.1, SNR G8. 31-0.09, and Suzaku Src2. 

SNR G8.7— 0.1: We discuss the association between 
HESS J1804-216 and SNR G8. 7-0.1. Figure [6] shows 
the leptonic model curves used to fit the TeV gamma-ray 
spectra obtained by H.E.S.S. and CANGAROO-III. We 
examined plausible inverse Compton (IC) model curves 
on the assumption that the electron spectrum is a single 
power-law with an exponential cutoff, 

dN, 
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Fig. 4.— Spectral energy distribution (SED) of HESS 
J1804-216, SNR G8.7-0.1, SNR G8.31-0.09, PWN of PSR 
B1800— 21, Suzaku Srcl, and Suzaku Src2 in all energy bands. 
The data points derived from this v^^ork are represented by fiUed 
squares and references to others are given in Tabled 
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where is the normalization factor, E„ is the electron 



Fig. 5. — SED in the X-ray band. The dashed and dotted 
arrows show the upper limits for the diffuse source obtained by 
Suzaku and Chandra, respectively. The dashed, dot-and-short- 
dashed, dot-and-long-dashed, and solid closed regions show the 
error regions of Suzaku Srcl, Suzaku Src2 with Suzaku telescope, 
Suzaku Src2 with Chandra telescope, and PWN of PSR B1800-21 
with Chandra telescope, respectively. 

energy, Fg is the spectral index of the injected electrons, 
and Emax_e is the maximum electron energy. In this 
paper, we do not consider hard spectra of electrons with 

Te < 1. 

To obtain the entire IC model curve, we used IR and 
optical (starlight) photon fields for the target photons 
in addition to the cosmic microwave background (CMB) 
field density. Here, we used an interstellar radiation field 
(ISRF) d erived from the lat est (v50p) GALPR OP pack- 
age (Port er fc Stron g 2005; Strong fc Moskalen ko 2006). 
The ISRF was given for three components (CMB, IR 
from dust, and optical starlight) as a function of the 
distance from the Galactic center, R (in kpc), and the 
distance from the Galactic plane, z (in kpc). We ex- 
tracted the spectra at (i?, z) = (3.8, -0.01) at the SNR 
G8.7— 0.1 position, as shown in Fig. [71 The X-ray upper 
limits for diffuse emission constrained Fg to be less than 
2.5 for the IC model curves. We obtained Emax_e to fit 
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the IC model curves to the TeV gamma-ray spectra for 
each fixed spectral index Fg (Table [2]). 

The electrons causing the IC scattering also emit syn- 
chrotron radiation. In the radio band, the spectrum of 
the whole region of G8.7— 0.1 is regarded as the upper 
limits. Here, the upper limits in the near-infrared (NIR) 
band, in Fig. 21 are neglected because they are for point 
sources. The magnetic field is constrained by the ra- 
dio upper limits and the X-ray upper limits for diffuse 
emission. Table [5] lists the upper limits of the magnetic 
field, maximum electron energy, and the total energy of 
electrons above 0.51 MeV, 

We = / dE, , (2) 

Jo. 51 MeV "-^e 

at d = 4.8 kpc for each electron spectral index Fg. We is 
estimated to be less than ~ 10^" ergs for Fg < 2.5. Al- 
though the total explosion energy of some supernovae are 
much higher, for example, the total energy of SN20031w, 
a type Ic supernova , is estimated to be ~ 6 x 10^^ ergs 
(|Mazzali et al.l [20061 ) . the typical total energy of a su- 
pernova explosion is estimated to be ~ 10^^ ergs. If the 
efficiency of the energy to accelerate the electrons is 10%, 
the obtained We can satisfy the SNR origin scenario. 

Additionally, the hadronic scenario can also explain 
the TeV gamma-ray spectrum for G8.7— 0.1. In Fig. [3 
we show the i:^ decay model curves with the assumption 
that the proton spectrum is a single power-law with an 
exponential cutoff, 

^ = KpEf- exp [-Ep/Er,,a._p) , (3) 

where Kp is the normalization factor, Ep is the pro- 
ton energy, Fp is the spectral index of the injected pro- 
tons, and Emax_p is the maximum proton energy. In 
this paper, we do not consider hard pr oton spectra with 
F„ < 1. The EGRET up per limit (|Kargaltsev et all 
l2007bt iHartman et al.l I1999D constrained Fp to be less 
than 2.3. We obtained Emax_p to fit the 7r° decay model 
curves to the TeV gamma-ray spectra for each fixed spec- 
tral index Fp. Table[3]lists the fitting parameters and the 
total energy of protons above 1 GeV, 

r AT 
Wp = / Ep-J^ dEp , (4) 

Jl GeV "-^P 

at a distance for SNR G8.7— 0.1 of 4.8 kpc with the as- 
sumption that the interstellar medium (ISM) density is 
n = 1 cm^'^. 

The total energy of protons, Wp, is 9.3 x 10^'^ {d/ 4.8 
kpc)^(n/l cm~'^)~^ergs with Fp = 2.3. For a typi- 
cal total supernova explosion energy of ~ 10^^ ergs, if 
the efficiency of the energy to accelerate the protons is 
10%, the ISM density should be n - lO^cm^^. With 
the assumption of Fp = 1.0, Wp is 9.6 x l0^^{d/4.8 
kpc)^(ri/l cm~'^)~^ ergs, and the ISM density should 
be n ~ 10 cm~'^ for an efficiency of 10%. However, 
iFinlev fc Ogelmanl () 19941 ) derived the electron density 
Ue to be (0.1 - 0.2)(d/6 kpc)-!/^ ^m-^ for an X-ray 
emitting gas in the remnant, based on the iJO^S^T" ob- 
servation. In addition, the lack of the X-ray emis- 
sion in the TeV gamma-ray emission region, shown in 
Fig. [2l indicates that the ambient density at the region 



is lower than that at the X-ray emitting region, that is, 
n < 0.1 — 0.2 cm"'^. If it were, the hadronic scenario 
would be unlikely On the other hand, a GMC, SRBY3, 
is near (in projection) to HESS J1804-216. The offset, 
the angular radius, the distance, and the mass of the 
cloud are ~ 0.1°, 0.2°, 5.3 kpc, 55 .6 x IO^Mq, respec- 
tively (jCramphorn fc Sunvaevll2002r ). The mean number 
density is calculated to be ~ 4.5 x lO^cm"'^ under the 
assumption that it is spherical symmetry. If it is associ- 
ated with HESS J1804— 216, there can be a high-density 
medium, and the hadronic scenario can be satisfied. 

Additionally, we estimated Kpe, the number ratio of 
protons to primary electrons. Under the assumption that 
the spectral indices of the protons and primary electrons 
are the same, Fp = Fg, Kpe is derived from eqs.([T])([3]) as 
Kpe — Kp/ Ke- Generally the maximum electron energy 
E„iax_e is the same as Emax.p, or Emax.e is lower than 
Emax_p due to a coohng effect, that is, E^ax.e < E^ax.p- 
In Fig. [51 we show the leptonic model curves with the 
assumption of Kpe ^ 10^(n/l cm~'^)~^, B = 3 jiG, Fg — 
2.0, and E^ax^e = Emax.p — 16 TeV. The leptonic model 
curves move up with the lower Kpe. With the assumption 
of B = 3/iG, Kpe is constrained to be Kpe > 10^(n/l 
cm^^)^^, because the IC model curve cannot exceed the 
TeV gamma-ray spectrum. If a GMC is associated with 
HESS J1804-216, and n ~ 10^ cm'^^ the obtained Kpe 
is consistent with that of the average cosmic rays in our 
Galaxy, Kpe ^ 10^. On the other hand, for a high KpeU 
of ~ 10"*, we need to consider the contribution of the 
emissi on from secondary electrons produced by charge d 
pions ([Yamazaki et al.ll2006t [Pfrommer &: EnBlinll2004l ). 
However, a detailed discussion about secondary electrons 
is beyond the scope of this paper. 

Suzaku Src2: For Suzaku Src2, we used the overlap- 
ping error region obtained by Suzaku and Chandra. In 
Fig. [9l we examined plausible synchrotron model curve 
within the error of the unabsorbed X-ray flux on the 
same assumption as eq. (jl} with adopting B—3 /iG and 
Fg — 2.0. The obtained IC model curves for the same 
electron spectrum cannot reproduce the TeV gamma-ray 
spectrum. Even if the higher IR energy density than that 
in GALPROP is adopted, the cutoff energy of the mod- 
els does not change from ^10 TeV, while the observed 
cutoff energy should be < 0.3 TeV; such a low cutoff en- 
ergy is due to the soft spectrum of HESS J1804-216. A 
stronger magnetic field can reduce the cutoff energy of 
the IC model curves. In order to make the cutoff energy 
of 0.3 TeV for Suzaku Src2, the magnetic field should be 
more than 7 mG, and to fit the TeV gamma-ray spectrum 
with the magnetic field the IR energy density should be 
-10^ times higher than that in GALPROP. These val- 
ues are unlikely. We therefore conclude that the TeV 
gamma-ray spectra are not produced by the IC model 
curves for the electron spectrum. 

We also examined the bremsstrahlung emission for the 
same electron spectrum. To fit the TeV gamma-ray spec- 
trum, the ISM density should be more than ~ 10® cm~^, 
and the magnetic field should be stronger than 1.8 mG 
for Suzaku Src2. Due to such incredibly high values of 
the ISM density and the magnetic field, we also reject 
the scenario that the bremsstrahlung radiation produces 
the TeV gamma-ray spectrum. 
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Fig. 6.— Leptonic scenario for SNR G8. 7-0.1. The thick soUd 
curves show the synchrotron emissions, and the thin soUd curves 
show the total IC spectra with fixed spectral indices of Fe = 1.0 (a), 
Fe = 2.0 (b), and Fe = 2.5 (c), and with B of 50^*0, 8/iG, 3/iG, 
respectively. The dotted, dot-dashed, and dashed curves show the 
IC spectra (a) on each of CMB, IR and starlight, shown in Fig.[7l 
respectively. The dash-dot-dotted curve shows the 1 year, 5(t sen- 
sitivity for the GLAST LAT taking into acc ount the d iffuse back- 
ground at the position of HESS J1804-216 \GL AST h AT: 2007). 

A simple solution to explain the observed TeV gamma- 
ray spectrum is that accelerated protons produce the 
gamma rays, as shown in Fig. [9l Therefore, the TeV 
gamma rays and the X-rays could be produced by ac- 
celerated protons and electrons, respectively. Assuming 
the spectral index, Tp = = 2.0, the maximum elec- 



tron energy is E„ 



max^p 



16 TeV (Table O, 



and this maximum energy makes the lower limit of the 
magnetic field oi B > 130 /zG to fit the X-ray data. 
With the magnetic field, we calculated the lower limit 
of Kpe > 6.6 X 10''(n/l cm^'^)^^. This value is much 
higher than that of the average cosmic rays in our Galaxy, 



pe 



10^ 



G8.31-0.09: We consider the possibility that HESS 
J1804-216 is associated with G8.31-0.09. For the lep- 
tonic scenario, it is difficult to explain the fact that the 
size of G8. 31— 0.09 in the radio band is much smaller than 
that in the TeV emission region. On the other hand, if 
the accelerated protons produced the TeV gamma rays, 
the difference in the size could be explained due to the 
difference of the diffusion length between protons and 
primary electrons. However, we need additional informa- 
tion, either on the distance or from other wavelengths, 
to discuss the energetics or Kpe- 

5.2. PWN origin 

At the present time, PWNe are the largest population 
of identified TeV galactic sources, and the number of TeV 
gamma-ray sources that are lo cated near a p lausible pul- 
sar candidate are increasing ()Hintonll200l . Therefore, 
it is very important to consider the PWN origin. The 
PWN candidates associated with HESS J1804-216 are 
the PWN of PSR B18G0-21, and Suzaku Src2. 

PWN of PSR B1800-21: To estimate the total en- 
ergy supplied from a pulsar, we should consider any 
braking effects for spin-down luminosity. If a pulsar 
spins down from an initial spin period of Pq according 
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Fig. 7. — Interstellar radiation field from the GALPROP pack- 
age (v50p) at {R,z) = (3.8,-0.01) kpe for G8.7-0.1. From lower 
frequencies, CMB (dotted), IR emission from interstellar dust (dot- 
dashed), and optical photons from stars (dashed) are shown. The 
solid line represents the sum of the three components. 
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Fig. 8. — Hadronic scenario for SNR G8. 7-0.1. The thin solid 
curves show the tt" decay model curves with fixed spectral indices 
of Fp = 2.3, Fp = 2.0, and Fp = 1.0, respectively. The thick soUd 
curve shows the total IC spectra, and the dot-and-long dashed 
curve shows the synchrotron emissions from the primary electrons 
with the assumption of Fe = 2.0, B = 3fiG, and Kpe = 10^(n/l 
cm~'^)~^. The dotted, dot-dashed, and dashed curves show the 
IC spectra on each of CMB, IR and starlight, shown in Fig. [7] 
respectively. 

to 17 = —kfl"'"' assuming both and k are constant, 
where ntr is the braking index of the pulsar, the spin- 
down luminosity, E(t), is given as 



To 



ribr - 1 



(5) 



To 



Pit) 



(ntr ~ l)Po {ntr - l)P(i) 



t . 



(6) 



where tq is the initial spin-down timescale 
([Pacini &: Salvatil Il973l ). Therefore, the total energy 
that the pulsar has lost over its age, r, is 



B — 



1 — a 
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Fig. 9. — SED with Icptonic model curves to fit Suzaku Src2 
and tt" decay model curves to fit TeV gamma rays. The thick 
solid curve shows the total IC spectra, and the dot-and-long dashed 
curve shows the synchrotron emissions to fit Src2 with the assump- 
tion of a spectral index of Fe = 2.0 and B = 3 fiG. The dotted, dot- 
dashed, and dashed curves show the IC spectra on each of CMB, 
IR and starlight, shown in Fig. [T] respectively. The long dashed 
curve shows the bremsstrahlung for the same electron spectrum 
with the assumption of n = l.Ocm"'^. The thin solid curves show 
the tt" decay model curves with fixed spectral indices of Fp = 2.3, 
Fp = 2.0, and Fp = 1.0, respectively. 
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Fig. 10. — SED and the leptonic model curves for the PWN of 
PSR B1800— 21 with a time-dependent rate of electrons-injection 
with T = 16 kyears, tq = 700 years, a braking index of nj^ ^ 1.6, 
Fe = 1.5, and B = 8 fj,G. The right-hand curves show the IC com- 
ponent, and the left-hand curves show the synchrotron emissions. 
The dot-dashed curves show the spectra produced by old electrons 
(t < 1000 years), and the dotted curves show the spectra produced 
by young electrons (t > 15000 years). The thin solid curves repre- 
sent their medium per 1000 years (1000 < t < 15000 years), and 
the thick solid curves show their total. 

For the Crab pulsar, produced in the supernova of 1054, 
a To of 700 years is derived from eq. ([6]), and the total 
energy is derived as i?pwn = 3.8 x 10^^ ergs with a brakin g 
index of Ubr — 2.5 from eq. ([7]) (Livi ngstone et ani2006D . 
For PSR B1800— 21, the braking index and the true age 
are unknown. The measured braking i ndices of other pul- 
sars f all in the range of 1.4 < Ubr < 3 (jLivingstone et al] 
|2006( ). Therefore, we estimated the total energies for sev- 
eral braking indices with E{t) — 2.2 x 10'^^ ergs/s, an age 
of r = 15.8 kyears (the same value as the characteristic 
age), and tq = 700 years (the same as the Crab pulsar), 




Energy [eV] 



Fig. 11. — SED and the leptonic model curves for Suzaku 
Src2 with a time-dependent rate of electrons-injection with r = 16 
kyears, tq = 700 years, a braking index of n^^ ~ 1.8, Fe = 1.5, and 
B = 8 fiG. The model curves are the same as in Fig. 1101 
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Fig. 12.— Hadronic scenario for PWN of PSR B1800-21 (a) 
and Suzaku Src2 (b). The thin solid curves show the tt" decay 
model curves with fixed spectral indices of Fp = 2.3, Fp = 2.0, and 
Fp = 1.5, respectively. The thick solid curves show the total IC 
spectra, and the dashed curves show the total synchrotron emis- 
sions from primary electrons with r = 16 kyears, tq = 700 years, a 
braking index of n^r ~ 1.6, Fe = 1.5, B = 4mG for PWN of PSR 
B1800-21 (a), and B = 600^0 for Src2 (b). 

which are listed in Table S) 

For the leptonic model, we assumed that the rate of 
electrons- injection, K'^{t), proportionally decreases with 

E{t), that is, K'^(t) oc E{t), and we considered the 
cooling effect due to synchrotron and IC energy losses 
(|Funk et al.ll2007[ ). The old electrons that were injected 
with high-power initial spin-down luminosity have been 
cooled. On the other hand the young electrons have 
less time to be cooled. Figure [TU] shows the SED with 
the leptonic model curves with a time-dependent rate of 
electrons-injection with tq — 700 years, a braking index 
of Ubr ~ 1-6, an age of r = 16 kyears, Fe = 1.5, and 
B =8 /iC: To obtain the model curves, we used the bro- 
ken energy due to synchrotron cooling and IC cooling, 
Ebr{t) ~ 3mlcyi4aTiT-t){By{8TT) + Wph)), where ar 
is the Thomson cross-section and Wph is the total en- 
ergy density of the photon field. For the PWN of PSR 
B1800— 21, we used the error region for the entire PWN 
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(inner + outer) obtained by Chandra (jKargaltsev et all 
I2007af ). The old electrons' are responsible for the TeV 
gamma-ray emission, and the young electrons are respon- 
sible for the X-ray emission. This picture may also ex- 
plain the size difference between the TeV gamma-ray and 
the X-ray emission regions, because the old electrons can 
extend further than the young ones. To calculate the to- 
tal energy of electrons. We, for the PWN origin scenario, 
we substituted Ei,r{t) for E^ax^e- Then, We is estimated 
to be ~ 1.7x 10^^ ergs for the above parameters. An E'pwn 
of 1.3 X 10^^ ergs for the same parameters can easily ex- 
plain the amount of electron energy, and the efficiency of 
the energy to accelerate electrons is ~ 1.3 x 10~^%. Note 
that the above parameters are not unique, and other 
combinations of a tq and a braking index are possible. 
If we adopt tq = 30 years, acceptable model curves can 
be derived with a braking index of nt,r ^ 2.2, an age of 
T = 16 kyears, Ve = 1-5, and B = 8 iiG. For the lat- 
ter parameters, E'pwn is estimated to be 2.3 x 10^^ ergs, 
and We is ~ 1.7 x 10"** ergs. The efficiency of the energy 
to accelerate the electrons is ^ 7.4 x 10~^%. In either 
case, a quite high value for the total energy, -Epwii, that 
the pulsar has lost is required to explain the observed 
spectra. 

We consider the hadronic model given in Fig. [T^ 
Table |3] lists the fitting parameters and the total en- 
ergy of protons Wp at a distance for PSR B1800— 21 
of 3.84 kpc with the assumption that the ISM density 
is n = 1 cm~^. Wp with Tp of 2.3 is estimated to be 
5.9 X 10^\d/3M kpc)2(n/l cm-3)-iergs. With a Tp of 
1.0, Wp is estimated to be 6.1 x I0^°{d/3M kpcf{n/l 
cm~^)~^ergs. The magnetic field is constrained to be 
B > 4mG with the assumption of Fp = 1.5 and the con- 
straint of Ernax_e < E^ax.p = 7.9 TcV (Table [3]). The 
magnetic field is very high. Kpe is estimated to be more 
than 1.6 x 10^(n/lcm-3)-^ If a GMC is associated with 
HESS J1804-216, and n - 10^ cm-^, the obtained Kpe 
is consistent with that of the average cosmic rays in our 
Galaxy, Kpe ^ 10^. However, the lower limit of the mag- 
netic field is much higher than the typical interstellar 
magnetic field, and there is no evidence supporting such 
a high magnetic field in this position. 

Suzaku Src2: For the Icptonic model, a similar sce- 
nario as that for PWN of PSR B1800-21 can be con- 
sidered for Suzaku Src2. However, because we have no 
information about the age and the spin-down luminosity 
of Suzaku Src2, we assume the same age and spin-down 
luminosity as that of PSR B1800-21. Figure fTTl shows 
leptonic model curves for Suzaku Src2 with tq = 700 
years, a braking index of n^r ~ 1-8, an age of r = 16 
kyears, Ve — 1-5, and B = 8 jiG. The total energy of 
the electrons is calculated to be We = 1.8 x 10''^(d/3.84 
kpc)^ ergs, and E'pwn is 1-2 x 10^^ ergs. This E'pwn would 
be proportionally reduced if E{t) is lower than that of 
PSR B1800— 21. For example, E'pwn is estimated to be 
~ 10^" ergs with the assumption of E{t) ^ 10^^ ergss"^. 
Radio and X-ray observations with higher time resolu- 
tions might be able to detect a sub-second pulse period 
from this source, and further observations in the IR to 
UV band for a diffuse source in the extended region of 
HESS J1804— 216 could contribute to a tighter solution 



for the leptonic models. 

We consider the hadronic model as given in Fig. [T2l 
The magnetic field is constrained to be i? > 600 jiG 
with the assumption of Fp = 1.5 and the constraint of 
Emax_e < Emax_p = 7.9 TcV (Tablc[3l). Kpe is estimated 
to be more than 1.5 x 10'^(n/lcm^'^)^^. The lower limit 
of the magnetic field is much higher than the typical in- 
terstellar magnetic field. 

6. CONCLUSIONS 

CANGAROO-III observed HESS J1804-216, and de- 
tected gamma rays above 600 GeV at the lOcr level during 
an effective exposure of 76 hr. The obtained differen- 
tial flux is consistent with the previous H.E.S.S. result, 
and the obtained morphology shows extended emission 
compared to our Point Spread Function. We have dis- 
cussed the radiation mechanism and considered the pro- 
posed counterparts to the TeV gamma-ray source. For 
the SNR scenario, the most plausible counterpart is the 
SNR G8.7— 0.1, and both hadronic and leptonic processes 
can produce the observed TeV gamma-ray spectrum. We 
obtained upper limits of spectral indices, Fg < 2.5 for 
leptonic scenario, and Fp < 2.3 for hadronic scenario. 
The total energy of the primary electrons, We, is esti- 
mated to be 3.3 X 10^"(d/4.8 kpc)^ ergs with F^ = 2.5. 
The total energy of the protons, Wp, is estimated to be 
9.3 X 10'5i(d/4.8 kpc)2(7i/l cm-3)-iergs with Fp = 2.3. 
The obtained Wp indicates that a molecular cloud might 
be associated to satisfy the hadronic scenario for the en- 
ergetics. Extending this to include an estimation of the 
emission from secondary electrons will provide a tighter 
solution for hadronic scenario. For the PWN scenario, we 
discussed both leptonic and hadronic processes. We con- 
firmed that the leptonic model with a time-dependent 
rate of electron-injection while considering the braking 
effect for the spin-down luminosity and the cooling ef- 
fect due to synchrotron and IC energy losses could ex- 
plain both the high TeV gamma-ray flux and the low 
X-ray flux. However, a quite high value for the to- 
tal energy that the pulsar has lost is required for PSR 
B1800-21 in this model, Ep^n ~ 10^^ ergs. For Suzaku 
Src2, we could obtain an acceptable model curve, and 
^-pwii was estimated to be ~ 10^° ergs with the assump- 
tion of E{t) ~ 10'^^ ergss"^. However, we need radio or 
X-rays observations to determine whether or not there is 
a pulse period in order to discuss the energetics further. 

The model curves that are given in this paper indicate 
that GLAST could determine the spectral index of the 
accelerated particles, and further observations in the IR 
to UV band for a diffuse source could give a solution of 
the counterpart and the radiation mechanism. 

The authors would like to thank B. Lott for providing 
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ence and Technology (MEXT), the Australian Research 
Council, and the Inter-University Research Program of 
the Institute for Cosmic Ray Research. The work is also 
supported by a Grant-in- Aid for the 21st century center 
of excellence programs "Center for Diversity and Uni- 
versality in Physics" and "Quantum Extreme Systems 
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TABLE 1 

Summary of data used in the SED analysis. 



Instrument 


Source name Symbol 


Kcierence 


VLA &nd other r&dio 


G8. 7—0.1 filled circle 




VLA 


G8. 31— 0.09 open circle 


(2) 


Suzaku XIS 


Suzaku Srcl closed region (dashed) 


(3) 


Suzaku XIS 


Suzaku Src2 closed region (dot-and-short dashed) 


(3) 


Suzaku XIS 


diffuse UL (90% C.L.) arrow (dashed) 


(3) 


Chandra ACTS 


Suzaku Src2 closed region (dot-and-long dashed) 


(4) 


Chandra ACIS 


PWN of PSR B1800-21 closed region (soUd) 


(5) 


Chandra ACIS 


diffuse UL arrow (dotted) 


(4) 


INTEGRAL IBIS 


Upper limit arrow (cross) 


(4) 


EGRET 


Upper limit arrow (solid) 


(4)(6) 


GLAST 


1 year sensitivity dash-dot-dotted line 


(7) 


2MASS 


NIR point UL triangle 


(8) 


H.E.S.S. 


HESS J1804-216 open square 


(9) 


CANGAROO-III 


HESS J1804-216 filled square 


this work 



References. — (1) Kassim & Weiler 1990b; (2) Brogan et al. 2006; (3) Bamba et al. 2007; 
(4) Kargaltsev et al. 2007b; (5) Kargaltsev et al. 2007a; (6) Hartman et al. 1999; 
(7) GLAST LAT 2007; (8) Skrutskie et al. 2006; (9) Aharonian et al. 2006a 



TABLE 2 

Summary of parameters used in the 

LEPTONIC model SHOWN IN FiG. [6] 





2.5 


2.0 


1.0 


Emax.e [TcV] 


3.6 


2.3 


1.3 


We [lO-lSergs] 


330 


17 


1.4 


Upper limit of B [ iiG ] 


1 


8 


50 



TABLE 3 

Summary of parameters used in the 
hadronic model. 



TABLE 4 
Braking index dependence of the 
parameters. 



braking index 


-Epwn [orgs] 


Po [ms] 


3 


2.6 X 10'*'' 


210 


2.5 


5.8 X 10*9 


120 


2.2 


1.3 X lO^'' 


72 


1.8 


1.2 X 10^1 


19 


1.6 


1.3 X 10^2 


5.2 


1.4 


1.7 X lO^-i 


0.37 





2.3 


2.0 


1.5 


1.0 


Emax^p [TeV] 


30 


16 


7.9 


5.4 


Wp" [10^1 ergs] 


9.3 


3.1 


1.3 


0.96 


VKp''[105iergs] 


5.9 


2.0 


0.84 


0.61 



For SNR G8.7-0.1, d = 4.8 kpc. 
For PSR B1800-21, d = 3.84 kpc. 
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